The Phase Diagram of the Quark-Meson 

Model 

B.-J. Schaefer a and J. Wambach a ' & 

a Institut fur Kernphysik, TU Darmstadt, D-64289 Darmstadt, Germany 
b Gesellschaft fiir Schwerionenforschung GSI, D-64291 Darmstadt, Germany 

February 4, 2008 



PACS: 



Abstract 



Within the proper-time renormalization group approach, the chi- 
ral phase diagram of a two-flavor quark-meson model is studied. In 
the chiral limit, the location of the tricritical point which is linked 
to a Gaussian fixed point, is determined. For quark chemical poten- 
tials smaller than the tricritical one the second-order phase transition 
belongs to the 0(4) universality class. For temperatures below the tri- 
critical one we find initially a weak first-order phase transition which 
is commonly seen in model studies and also in recent lattice simula- 
tions. In addition, below temperatures of T < 17 MeV we find two 
phase transitions. The chiral restoration transition is initially also of 
first-order but turns into a second-order transition again. This leads 
to the possibility that there may be a "second tricritical" point in the 
QCD phase diagram in the chiral limit. 
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1 Introduction 



Recently, the phase diagram of the chiral phase transition has received re- 
newed interest with the emergence of improved lattice techniques at finite 
quark chemical potential, \i. Through heavy-ion collision experiments at 
RHIC, the LHC and the future GSI facility parts of the phase diagram may be 
accessible. As one of the prominent features of the QCD phase diagram, the 
existence of a (tri-) critical endpoint of a first-order phase transition line has 
emerged. This is strongly suggested by model calculations and also indicated 
by lattice simulations. The precise location of this point is still unknown. In 
the chiral limit, the critical endpoint becomes a tricritical one corresponding 
to a Gaussian fixed point. In this limit, the chiral finite-temperature transi- 
tion at vanishing \x is likely to be of second-order for two-flavor QCD, where 
the critical transition line is expected to fall into the 0(4)-universality class 
in three dimensions [L. For finite current masses, the 0(4)-transition line 
converts to a smooth crossover and the tricritical point turns into a fixed 
point belonging to the 3d Ising universality class [21 EH E] • 

Second-order phase transitions are characterized by long-wavelength fluc- 
tuations of the order parameter. Thus, the critical and non-critical regions 
of the phase diagram, including first-order phase transitions, can best be de- 
scribed by renormalization group (RG) methods. We employ an analytical 
RG approach, the proper-time renormalization group (PTRG) El U\ > to 
investigate the QCD phase diagram for two quark flavors. As a result of 
integrating the quantum and thermal fluctuations with a certain cutoff scale 
k, we obtain a flow equation for the scale-dependent action. Following the 
scale evolution towards k = 0, we arrive at the effective action which contains 
all thermodynamically relevant universal and non-universal information. 

The paper is organized as follow: In the Sec. 121 we derive the analyti- 
cal PTRG flow equation for an optimized smearing function. This allows 
to find analytical threshold functions for finite temperatures and chemical 
potentials. In this way, the numerical summation of Matsubara terms can 
be circumvented and the numerical error especially for small temperatures 
during the evolution is reduced. In Sec. 12.11 we consider analytical limits of 
the full flow equation in order to elucidate its intuitive structure and sketch 
the numerical implementation in Sec. 12.21 The resulting phase diagram for 
massless quarks is presented and explored in some detail in Sec. EJ Finally, 
we summarize and conclude in Sec. 
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2 The flow equations for finite T and \i 



As an effective realization of two-flavor QCD, we use the effective action of a 
linear quark-meson model defined at a scale A in the ultraviolet (UV) region 
of the theory in four Euclidean dimensions by 1 

r A [$] = jd A x jg[7 M fy+0(*+"*7?75)]g + \{d^af + ^vf) 2 + V A (a 2 + n 2 



(1) 

with one isoscalar-scalar a-meson and three isovector-pseudoscalar pions 7?. 
The fermionic fields q(x) and q(x) incorporate Nf = 2 flavor and N c = 3 
color degrees of freedom, f labels the three Pauli matrices. The Yukawa 
coupling g is kept constant during the RG evolution. The purely mesonic 
general effective potential term V A is a function of the 0(4)-symmetric field 
<f) 2 = a 2 +n 2 . As initial condition in the UV we choose for V A the truncation 

V A (a 2 + if 2 ) = ^(a 2 + n 2 ) + ^V + 7 ? 2 ) 2 (2) 

with the parameters m\ and Aa to be specified later. The wavefunction 
renormalizations for the meson- and quark fields are neglected in this ap- 
proximation. 

The coarse-grained, infrared-finite and renormalizat ion-group improved 
flow equation for the scale dependent effective action r fc [$] of a given theory 
is governed by 

00 

d t T k [<S>] = ~\\% [dth(rk 2 )} Trexp (-rif } [$]) , (3) 


where Tr denotes a four- dimensional momentum integration and a trace over 
all given inner spaces (e.g. Dirac, color and/or flavor-space). Details con- 
cerning the derivation of this flow equation and inherent approximations can 
be found in Refs. j3J |H] . 1^ [$] represents the full inverse propagator and is 
given by a second functional derivative of with respect to the appropriate 
field components 

rf [•] = (4) 



1 The argument $ of the effective action is a short- hand notation of all incorporated 
fields. 
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If one replaces the effective action derivative r^ 2 ^ on the right-hand side 
of Eq. (J3J) by the corresponding classical action one obtains the one-loop 
effective action 0[7]. The replacement of the classical action by the 
scale-dependent effective action on the right-hand side of Eq. Q is called 
RG improvement. 

By means of a Schwinger proper-time regularization, an infrared-finite 
renormalization group equation (PTRG) is obtained. The infrared (IR) scale 
is given by t = ln(fc/A) where A denotes the high momentum (UV) scale. 
One important ingredient of the flow equation Q is the a priori unknown 
smearing or blocking function /^(rfc 2 ) which regulates the UV and IR diver- 
gent proper-time integration, respectively. 

The convergence of the approximate solution of the flow equation is partly 
controlled by the regulator and can be accelerated if a certain blocking func- 
tion is chosen. In previous works [HJ EI| the regulator was given by 



for d = 4 dimensions. In Refs. [HJ EH details concerning the choice of regu- 
lators and the scheme dependence of the proper-time flow equations for the 
quark-meson model can be found. 

This choice has the disadvantage that the generalization to finite tem- 
perature results in Matsubara sums which have to be evaluated numerically. 
For higher smearing functions i.e. for functions with indices i > 1 more and 
more additional terms in the flow equations are generated at finite T and 
[i and render the flow equations more complicated. The choice of the same 
regulator but now for i = 1 and in d = 3 dimensions, produces an analytic 
RG flow equation in T as well as in /x and is given explicitly by 



Then only one threshold function for each physical degree of freedom arises 
and the structure of the full RG flow equation for the thermodynamic poten- 
tial in the vacuum becomes simple and physically intuitive. Also the finite 
T and \x generalization is straightforward. 

In thermal equilibrium, we apply the Matsubara technique and follow 
the standard procedure. The quark chemical potential is introduced in the 
fermionic piece of the Lagrangian by the replacement do — > do — i\i in the 
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(6) 
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time-component of the derivative. The finite-temperature extension of this 
PTRG method is described in detail in Ref. [H] for an 0(7V)-model and in 
Ref. JU] for the quark-meson model. 

In the Matsubara formalism, the integration over the zero- component of 
the momentum is replaced by a summation over Matsubara frequencies, u n , 
for mesons and, u n , for the quarks: 



2nnT 



Vn 



(2n + 1)ttT 



such that for an arbitrary function g(p) 



d 4 p 



d 3 p 



g(p, 



n G 



(7) 



(8) 



for Fermions and Bosons, respectively. 

In order to solve the RG equation (J3J), a truncation of the effective ac- 
tion is necessary. We will use the lowest-order in a derivative expansion of 
the effective action given in Eq. JIJ at a given UV scale A. This yields a 
single partial differential equation for the scale-dependent grand canonical 
thermodynamic potential Q(T, \i\ 0) which depends on T, /i and the fields </>: 



T 
2 



dr 

T 



dtf, 



(t=l,d=3) 
k 



d 3 q 



Tr 



(2tt) 3 

e -r[(^+iAi) 2 +g a +s 2 ^ 2 ] 



(9) 



Only the fermionic potential contribution is modified explicitly by the finite 
quark chemical potential. At this stage, one sees nicely the additive separa- 
tion of the mesonic and fermionic flow in this approximation. Note, that we 
use the same blocking function for the bosonic and fermionic flow-equation 
contributions. 

Evaluating the traces and performing the remaining three-momentum 
integrals, we finally obtain 



d t tt k (T,fx; 



T 



2(4tt) 3 / 2 J r 5 /2 



dr r 



dtf, 



(i=l,d=3) 



(10) 



4N c N f e~ T 
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where the ^-dependence of the potential derivatives on the right-hand side 
has been suppressed. The primed potential denotes the 2 -derivative, i.e. 
Q' k := dQk/d(p 2 etc. In the following, we define the pion-, a-meson and 
quark energies as 



E w = v^ 2 + 2Q' k , E a = y/k 2 + 2Q' k + , E q = ^/k 2 + g 1 ^ . (11) 



Evaluating the potential at the global minimum, <j) 0l the scale-dependent 
effective meson masses are defined as m 2 fc = {2Q! k +A(f) 2 Q! k )\ ( f )=( j )0 = 40q^(0 o ) 
and m 2 k = 2f^|0 =( £ o = 0. The dynamically generated (constituent) quark 
masses m q ^ = g4>o are proportional to the minimum since we do not consider 
the running of the Yukawa coupling constant. 

Finally, the PTRG flow equation for the scale-dependent grand canonical 
potential can be integrated analytically, resulting in 



This flow equation is the full RG equation in the approximation considered for 
finite temperature and chemical potential. Note, that since we have neglected 
the running of the Yukawa coupling, the approximation does not correspond 
to the local potential approximation (LPA). It is possible to rewrite Eq. (|12j) 
by means of occupation numbers (cf. Ref. [IHII2])- 

The different degrees of freedom contribute in an additive way to the flow. 
One recognizes the three (degenerate) pion-, the sigma- and the quark/anti- 
quark threshold functions in the square brackets. The flow equation has an 
overall scale factor k , the correct dimension of the potential in d = 4 dimen- 
sions, which can be seen explicitly by rewriting all threshold functions in a 
dimensionless form. Due to the fermion loops, the fermionic contributions 
enter with a negative sign and have a degeneracy factor of (2s + l)xJV c x Nf 
with s = 1/2. Only in the quark/antiquark threshold functions the quark 
chemical potential enters with the appropriate sign as it should be. 




Eg+JJ, 

2T 



) 



(12) 
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2.1 Analytical limits 

It is instructive to investigate analytical limits of the full flow equation ()12j) . 
The flow in the vacuum is given by 



<9 t fi fc (O,O;0) = 



I2tt 2 



3 1 AN c N f 
E T E a E q 



(13) 



Due to spontaneous chiral symmetry breaking, the quark energies E q are 
always positive and contribute to the flow. In the mesonic threshold functions 
poles can emerge, depending on the shape of the potential. For instance, the 
pion pole for <fi = is determined exactly by the condition Q' k = —k 2 /2. 
However, for an appropriate choice of initial conditions at the UV scale A, 
the poles are never reached during the /c-evolution. 

For finite T but vanishing /i, the difference between fermionic and bosonic 
Matsubara sums in d t flk{T, 0) vanishes in the low-temperature limit and we 
regain the vacuum flow analytically. For finite fi and T = 0, on the other 
hand, the flow is governed by the non-analytic equation 



k 



12tt 2 



3 1 4N c N f 
E w E a E q 



{1-9^-Eq)} 



(14) 



reflecting the existence of the Fermi surface. Thus, in the zero-temperature 
limit, the tanh-expressions in the fermionic threshold functions of Eq. (|12j) 
degenerate to 29{E q — fi). It is also possible to analytically perform the zero 
chemical-potential limit of Eq. (fT2*j) which again yields the correct vacuum 
flow since the fermionic flow contribution in the curly brackets tends to two. 

For hadronic matter in thermal equilibrium, two new energy scales T 
and n enter the evolution equation. The quark chemical potential influences 
the bosonic part d t Q^ oson only implicitly through the meson masses. If 
the bosonic fluctuations are nee: lected by setting <9^ oson = 0, one obtains 
(for a constant Yukawa coupling) standard mean-field theory results for the 
thermodynamic potential. The remaining quark contribution can then easily 
be integrated [13]. This feature is also seen in other RG approaches (e.g. in 
Ref. |I1]). As is well known, mean- field theory does not always give a reliable 
description of the phase transition near criticality as it neglects the role of the 
important (thermodynamic) fluctuations. We will come back to this point 
below. 

Rewriting the step-function as in Eq. (JHJ), one recognizes easily that the 
quark flow contribution splits into a vacuum and finite-chemical potential 
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part with opposite sign. For ji < E q = g4>o,k=0: the finite-/i part vanishes 
and there is no distinction to the vacuum flow. In the vacuum, spontaneous 
chiral symmetry breaking is mainly driven by the quark fluctuations, while 
meson fluctuations tend to restore chiral symmetry due to the relative sign 
in the flow Eq. (jl2j) . The relative sign between the vacuum and finite density 
part in the quark sector of Eq. (|14|) suppresses the quark fluctuations which 
are important for chiral symmetry restoration at high densities. 

2.2 Solving the flow equation 

In order to solve the flow Eq. (|12j) numerically one has in principle two 
possibilities: either one can discretize the unknown potential on a 2 - 
grid or expand the potential in powers of <p 2 around its minimum <pQ. The 
advantage of a potential expansion is that only a few coupled flow equations 
have to be solved, depending on the number of couplings in the potential. 
But for each higher order of the potential expansion, a new coupled beta 
function arises which increases the numerical effort drastically. A further 
disadvantage is that the potential, after the evolution, is only known at the 
(local) minimum (J)q. 

The big advantage of the grid solution is that one knows the potential 
not only for the minimum but also for arbitrary 2 . This is of importance 
for a first-order phase transition where two degenerate minima of the po- 
tential emerge. In order to describe the first-order transition correctly, the 
knowledge of all local minima is required. This is cumbersome in a potential 
expansion, except for some simple potentials. For finite pion masses, i.e. with 
an additional explicit symmetry breaking term in the potential, every mini- 
mum has always a finite value and the symmetry is never exactly restored. 
Also in this case, a precise determination of the critical temperature of a first- 
order transition is very difficult within an expansion scheme around only one 
minimum. The solution on a grid is, however, computationally very demand- 
ing since one encounters highly-coupled and rather large matrix equations in 
order to obtain reliable accuracy. Nonetheless, for the reasons given above, 
we opted for the grid solution but will also compare the findings with results 
from a potential expansion. 

In order to solve the flow equation ()12|) . we discretize the field (f> 2 for a 
general potential term on a regular grid. Details of the algorithm which has 
been adapted to include fermions, can be found in Ref. [DJ and references 
therein. The task is then to solve the resulting closed coupled set of flow 
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equations for the potential by starting the evolution in the ultraviolet and 
integrate towards k = for each grid point. In this limit, all quantum- and 
thermal fluctuations are taken into account. It turns out that the numerical 
convergence is faster if one works with the dimensionful flow equation instead 
of the rescaled one. The initial conditions are chosen at the UV scale A in 
such a way that they match physical vacuum quantities obtained in the IR. 
Predictions for finite temperature and density are then possible without any 
further parameter adjustments. As initial condition we use the tree-level 
parameterization of the symmetric potential 

(15) 

A 

at the compositeness scale A = k xSB . For the results presented below we 
choose A = 500 MeV, A = 10 and 0q = 0. The number of grid points is 
typically chosen around 60 — 100 for a ^-interval between and 100. We 
have tested, that the numerical results do not change by varying the number 
of grid points. The precise value of </>q is not important, since it has little 
influence on the vacuum results in the infrared where chiral symmetry is 
spontaneously broken. The Yukawa coupling is fixed to g — 3.2 in order to 
reproduce a constituent quark mass of the order of 300 MeV. These choices of 
initial parameters yield a vacuum pion decay constant of f n ~ 87 MeV in the 
chiral limit which is consistent with values obtained from chiral perturbation 
theory ^5] and previous works [§]. For finite and realistic pion masses, this 
value is shifted to f n ~ 93 MeV. 

3 The phase diagram for Nf = 2 

Recently, significant progress has been made in studies of the QCD thermo- 
dynamics with chemical potential by numerical simulations of lattice gauge 
theory. Different methods of extrapolation to the finite-// region are in rea- 
sonable agreement [T7| [TBI CH] concerning the pseudocritical line T c (/z) out 
to chemical potential values of several 100 MeV [213 El 122] • This allows for 
a direct comparison of lattice simulations with our nonperturbative method. 
Some recent results of the leading non-trivial term of quadratic order in the 
Taylor expansion of the pseudocritical line are listed in Tab. ^ and compared 
to the one obtained with the PTRG method. 



^a(0 2 ) = ^(0 2 ) 2 -^0 2 
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rp d 2 T c 
c dfj, 2 




PTRG (full) 


-0.23 




PTRG (0 4 ) 


-0.258 




Lattice (Taylor) 


-0.14(6) 




Lattice (/!/) 


-0.101 




NJL (mean- field) 


-0.4 





Table 1: The curvature of the Nf = 2 transition line for small quark chemical 
potentials compared with different methods. 



The curvature of the RG without any truncation, labeled by PTRG (full), 
is smaller compared to a 4 truncation of the potential within the same 
RG approach, denoted by PTRG (0 4 ). At the expected tricritical point, a 
first-order phase transition emerges where higher monomials in the potential 
expansion will become necessary. One the other hand, we know that the 
truncation works well for the finite temperature and zero chemical potential 
phase transition. It turns out that the 4 truncation breaks down at finite 
quark chemical potentials around \x ~ 100 MeV [23] . 

Note that the RG results are obtained for zero quark masses in contrast 
to the lattice simulations (see e.g. table 1 in Ref. ^H] for recent lattice 
pion masses in physical units). On the lattice (e.g. Ref. [T7|). realistic light 
quark masses have not yet been achieved in a Taylor series estimate of the 
reweighting factor up to quadratic order in /z. The first Taylor coefficient 
corresponds to the curvature and is labeled by (Taylor) in the Table. The 
equivalent coefficient obtained with the imaginary chemical potential method 
in denoted by (/ij) in the Table [18} . However, the lattice results of Ref. [T7] 
suggest that any dependencies of the curvature on the quark masses are 
weak. Compared to the chiral limit we see a difference of a factor two. 
A mean-field calculation within a NJL-model in the chiral limit yields the 
largest curvature. However, the curvature depends strongly on the chosen 
parameter set for the NJL model [2"4*j . 

In Fig. ^ we show the chiral phase diagram for the two flavor linear 
quark- meson model in the chiral limit. At /x = 0, a second-order phase 
transition with a critical temperature of T c ~ 142 MeV is found in which 
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the spontaneously broken chiral symmetry is restored. Note that the critical 
temperature is not a universal quantity and depends on the fitted pion decay 
constant in the vacuum. The phase transition and fixed point structure is 
in the 0(4) universality class and was investigated within the same PTRG 
framework in detail in Ref. 



> 
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Figure 1: The phase diagram of the two-flavor quark-meson model in the 
chiral limit. The tricritical points are denoted by an asterisk. For fi below 
the upper tricritical point the transition is of second-order (dashed line) and 
belongs to the 0(4) universality class. For temperatures smaller than T < 17 
MeV the single first-order phase transition (solid line) splits into two phase 
transitions. The left transition line is always of first-order while the chiral 
restoration transition line around the splitting point is initially of first-order 
(solid line) and turns into a second-order (dashed line) phase transition for 
smaller temperatures. 

For realistic pion masses the second-order phase transition is washed out 
and becomes a smooth crossover. This results in a larger "critical" temper- 
ature of the order of 180 MeV |2*o] . 

For finite chemical potential, the second-order phase transition with 0(4) 
critical exponents persists up to a tricritical point. The location of this tri- 
critical point is found to be at T c ,tri — 52 MeV and n c ,tri — 251 MeV for the 
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initial conditions chosen. Within a nonperturbative expansion scheme, simi- 
lar results are found in Ref. (2H|. This point belongs to a trivial Gaussian fixed 
point with mean-field critical exponents as we have verified explicitly [2*7] . 

For finite (realistic) pion masses, model calculations and lattice simula- 
tions suggest a range of T ~ 100 — 180 MeV and \ib ~ 50 — 700 MeV for the 
location of the critical point of QCD [2 . The critical mean-field region of 
the endpoint, where the mean- field theory of phase transitions is still valid, 
is further expected to be small |3j. Nevertheless, the qualitative features of 
the phase diagram which we address, do not depend on the chosen initial 
conditions. Thus, the existence of this point, the shape of the transition 
lines and its universality class is a prediction within the underlying model. 

At higher chemical potentials and smaller temperatures, the phase tran- 
sition changes initially to a single first-order phase transition. For critical 
temperatures below T s ~ 17 MeV and around /x s ~ 263 MeV, however, we 
observe a splitting of the transition line and two phase transitions emerge. 
The left transition line (see magnified panel in Fig. [IJ represents a first-order 
transition down to the T = 0-axis. At this transition, the order parame- 
ter jumps not to zero but to a finite value (Fig. |2J). The chiral symmetry 
remains spontaneously broken and is only restored for higher chemical poten- 
tials which produce the second (right) transition line. For temperatures just 
below the splitting point, i.e. for T < 17 MeV, at the right transition line we 
initially find a second first-order transition where the order parameter again 
jumps to zero and chiral symmetry gets restored. But for smaller tempera- 
tures close to the /x-axis the order parameter tends smoothly to zero. Due to 
a finite grid spacing it is difficult for such small temperatures to determine 
the precise order of the transition. It seems, however, that it is again of 
second-order. The determination of the corresponding universality class is 
postponed to future work (.23]. Nevertheless, we infer that there must be a 
second tricritical point in the phase diagram. 

In Fig. |2I the order parameter, which corresponds to the expectation value 
( VEV) of the cr-field, is shown as function of the chemical potential for three 
different temperatures in the splitting area of the phase diagram (Fig. [TJ. 
Around the splitting point at T = 15 MeV, the order parameter depends 
slightly on the chemical potential \x before the first transition, thus washing 
out the edge at the gap. The //-dependence becomes weaker for decreasing 
temperatures. For T = 1 MeV there is almost no difference between the 
finite-density and vacuum evolution and the order parameter stays at the 
vacuum value <fi ~ 87 MeV. For T = 15 MeV two gaps in the order parameter 
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Figure 2: The order parameter versus chemical potential for three different 
temperatures. 



corresponding to two phase transitions can be seen. The magnitude of the 
first gap is almost constant and independent of the temperature. For smaller 
temperatures (T < 8 MeV), the VEV as a function of /x, smoothly decreases 
to zero after the first transition and stays zero for larger chemical potentials. 
This is a signal that the system undergoes a second-order transition again. 
Finally, chiral symmetry is completely restored. It is interesting to observe 
that all curves intersect roughly at the same critical chemical potential of the 
splitting point, fi s , where the magnitude of the VEV is about one half of the 
vacuum value. 

The exact location of the splitting point and the first-order gap of the 
order parameter depends on the chosen vacuum quark masses. For larger 
vacuum quark masses e.g. of the order of 380 MeV, corresponding to a larger 
Yukawa- coupling of 4.2, the location of the splitting point moves down to 
smaller temperatures and larger chemical potentials. The gap in the or- 
der parameter of the first-order transition increases with the quark masses. 
Thus, varying the vacuum quark masses within say 250 — 450 MeV, the 
area bounded by the transition lines is reduced. Yet, the order of the phase 
transition lines are not altered. 

It is also instructive to see both phase transitions in the thermodynamic 
potential explicitly. In Fig. 01 the fc-evolution below the splitting point is 
shown for T = 9 MeV and ji = 254 MeV for different (small) IR scales. The 



13 




20 40 60 80 100 

4 [MeV] 



Figure 3: The scale evolution of the thermodynamic potential towards the 
infrared is shown as a function of the field <j) below the splitting point for 
T = 6 MeV and \i = 254 MeV. (Lowest line k = 8 MeV,A/c = 1 MeV 
between each line). Since the chemical potential is larger than the critical 
one, both minima of the potential (one at ~ 65 MeV and the other one at 
~ 87 MeV) are no longer equal. 

potential exhibits two finite minima. 2 This situation is definitely different 
compared to that in Fig. HJ Here the scale evolution of a typical second- or 
first-order transition is shown in the left- and right panel, respectively. For a 
second-order phase transition we have only one (global) minimum. In the IR, 
the potential becomes flat (constant) below the minimum and is almost scale 
independent above . The scale evolution barely effects the minimum. For the 
chosen initial conditions, this situation is encountered for chemical potentials 
below 250 MeV and temperatures between 142 - 52 MeV. A typical first-order 
transition is shown in the right panel of Fig. linear the critical point. Here, 
a second minimum at the origin emerges and the global minimum jumps 
to zero for larger chemical potentials. In all cases, the potential tends to 
a convex function but the scale evolution of the nontrivial minimum is very 
weak. This justifies to terminate the numerical evolution at a finite but small 

2 For any finite IR scale (which is numerically always realized) the potential is not 
exactly convex. Only in the IR limit k = 0, it becomes convex as can be seen in the Fig.0] 

M 
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Figure 4: Typical scale evolution towards the infrared of a second-order phase 
transition (left panel) and first-order transition (right panel). (Lowest line 
k = 2 MeV, Ak = 1 MeV between each line. See text for further details). 

IR scale (typically < 1 — 10 MeV) [2E|. 

The pressure p(T,p) = — Q k=0 (T, /i, <p ) evaluated at the global mini- 
mum 0o is normalized to zero in the vacuum. For finite chemical potential 
and vanishing temperature, it stays strictly zero up to the transition point. 
Therefore, the quark number density, n q = dp/ dp, jumps to a finite value. 
Afterwards it increases monotonically. Quantitatively, the quark number 
density jumps to 0.1 fm -3 at /i ~ 242 MeV which corresponds to a normal- 
ized baryon density Pb/Po ~ 0.18. These values are obtained for a vacuum 
quark mass of ~ 275 MeV and increase significantly towards more realistic 
values if the vacuum quark mass is increased. These observations strongly 
suggest that the system undergoes a gas-liquid transition to bound quark 
matter, where the chiral symmetry remains spontaneously broken. Thus, 
the "triangular region" in Fig. ^ consists of a bound mixture of massive 
quarks/antiquarks, interacting massless pions and massive sigma mesons. 

4 Conclusions 

In order to provide a reliable identification of the phase structure of strongly 
interacting matter, one must invoke nonperturbative methods such as the 
Wilsonian formulation of the renormalization group. We have employed a 
self-consistent RG approach with a smooth proper-time regularization. As 
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an effective realization of QCD for two flavors, a linear quark-meson model 
has been used with an arbitrary 0(4)-symmetric mesonic potential term. 
To explore the phase diagram in the (T, /z)-plane we have extended previous 
studies to finite temperature and quark chemical potential. The resulting RG 
equations for the potential term have been solved numerically in the chiral 
limit. The first non-trivial coefficient of the transition line which belongs to 
the 0(4) universality class, has been determined and compared with recent 
lattice simulations. The curvature in the chiral limit is larger by a factor two 
as compared to lattice analyses. The latter, however, employ unrealistically 
large u — d current quark masses. 

The location of the tricritical point, governed by a Gaussian fixed point, 
has been determined. This structure of the phase diagram is similar to 
the results of other works, where the second-order transition line ends in a 
tricritical point and turns into a single phase transition of (weak) first-order 
for decreasing temperatures j23 I2H] ■ This first-order transition persists up 
to zero temperature and finite chemical potential. 

For fi > fi c we also find initially a first-order transition which splits into 
two phase transitions for temperatures below 17 MeV. The left transition line 
is always of first-order. The value for the splitting point depends on the initial 
condition for dynamically generated constituent quark mass in the vacuum. 
Apart from the dependencies on the initial conditions, further uncertainties 
arise from the truncation scheme of the flow equation which may affect the 
exact values of the critical temperature and chemical potential. It is difficult 
to obtain quantitative error estimates in this context. The main focus of this 
work is not a quantitative prediction of some non-universal numbers but the 
qualitative feature of the phase diagram of the two-flavor linear quark-meson 
model. 

For larger quark masses the splitting point moves towards smaller tem- 
peratures but all qualitative features of the transition lines persist. Along 
the chiral symmetry restoration line below the splitting point, the transition 
is initially also of first-order but then turns into a second-order transition. 
This leads to the suggestion that the QCD phase diagram may have a second 
"tricritical" point. 

The triangular region in the phase diagram which is a new feature of 
our calculation, is not related to spinodal curves or metastable phases. It 
is an effect driven by fluctuations. In a mean-field calculation, we only find 
one single first-order transition in which chiral symmetry is restored (see 
also 29] ). The full RG treatment, however, renders the order parameter and 
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hence the constituent quark mass finite beyond the first-order transition line. 
At the transition the quark-number density jumps from zero to a finite value 
and increases monotonically with increasing \i. We interprete this transition 
as the liquid-gas transition to a state of bound quark matter, consisting of 
interacting constituent quarks. 

There are several directions in which the present analysis should be ex- 
tended. A obvious one is to include explicit chiral symmetry breaking via 
finite current quark masses j2H|. This will turn the second-order transitions 
into smooth cross overs. From lattice studies it is well known that the strange 
quarks mass plays a decisive role in the location of the critical endpoint. For 
infinitely heavy strange quarks, in fact, it lies at /i = 0! Thus, it is manda- 
tory to extend our calculations to Nf = 3. While this is a computational 
challenge it will provide more realistic results. Work in this direction is in 
progress. 
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